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Abstract 

We discuss the generic phase diagrams of pure systems that remain fluid near zero temperature. 
We call this phase a quantum fluid. We argue that the signature of the transition is the change of 
sign of the chemical potential, being negative in the normal phase and becoming positive in the 
quantum fluid phase. We show that this change is characterized by a phenomenon that we call 
condensation, in which a macroscopic number of particles is in their own many-body ground state, 
a situation common to Fermi and Bose gases. We show that the ideal Bose-Einstein Condensation 
fits in this scenario, but that it also permits the occurrence of a situation that we may call "Fermi- 
Dirac Condensation" . In addition, we argue that this phenomenon is also behind the development 
of superfluid phases in real interacting fluids. However, only interacting systems may show the 
change from a thermal fluid to a quantum one as a true phase transition. As a corollary, this 
argument shows the necessity of the appearance of a "supersolid" phase. We also indicate how 
these ideas may be useful in the description of of experiments with ultracold alkali gases. 

PACS numbers: 67.25.D- Superfluid phase (^He), 67.30.H- Superfluid phase ('^He), 03.75.Hh Static properties 
of condensates; thermodynamical, statistical, and structural properties. 
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I. INTRODUCTION 



The resolution of the thermodynamics of interacting systems at very low temperatures 
starting from first principles is an aged endeavour that has even attracted the attention 
of great minds in the last and present Centuries, sparked early by the striking properties 
of Helium[lHl] and renewed with the spectacular experiments involving ultracold Alkaline 
gases[3H7]. Initiating with the crucial experimental and theoretical works by KapitzajT], 
Allen and Wisener[2], London [H], Tisza[9] and Landau [10], that set the tone for the fun- 
damental advances in the understanding of quantum macroscopic systems, and continuing 
with the most solid theoretical advances by Bogolubov[Tl] and Bardeen-Cooper-Schriefer 
(BCS) theories [T3] of interacting Bose and Fermi gases, a picture has emerged in which 
it is clear that at very low temperatures and even moderate pressures, certain fluids undergo 
a phase transition to a superfluid phase. This phase is envisioned as a fluid with two parts, 
one a "normal" fluid and another, a "superfluid" fraction that does not contribute to the 
entropy [01 ^0\. At zero temperature the system reaches its ground state but remains fluid, 
being superfluid all of it. The recent observation of quantized vortices in alkaline vapours 
of both Bose p^HTT] and Fermi [18] atoms, also confirms the superfluid nature of those gases 
and, hence, of the attained quantum macroscopic phase. 

Although a satisfactory understanding of the nature of the corresponding quantum states 
near and at zero temperature has been made, as mentioned above, and a lot is known 
regarding general thermodynamic aspects of these fluids, see Refs. [T9H23] for instance, we 
would to point out here an aspect that seems to have gone unnoticed and this is the role 
of the chemical potential as an indicator of the transition to a superfluid state, related to 
the main idea that the superfluid fraction of the fluid does not contribute to the entropy. 
Despite all efforts to portray the appearance of superfluidity as a condensation phenomenon 
a la Bose- Einstein Condensation (BEC), we shall point out a striking property of ideal Fermi 
gases that serves as an illustration of the main hypothesis here advanced, namely, that the 
condensation phenomenon is the appearance of the many-body ground state of the particles 
in the condensate, and not of the occupation of a single particle state. We insist right away 
that this hypothesis is independent of Bose or Fermi statistics, but certainly, ideal BEC is 
included by this identification. Of course, the quantitative details of behavior of different 
fluids may strongly depend on the statistics and other peculiarities of the atoms or molecules 
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in question, but no so their overall thermodynamic behavior leading to a superfluid phase. 

The present article is based on a thermodynamic approach rather than one based on 
many body quantum mechanics. As a matter of fact, we lack such a theory to fully back the 
hypothesis put forward here, and a hope is that these ideas may show the way for finding it. 
In the last sections we point out that the present experiments with trapped alkaline gases 
may be interpreted in such a way as to being a direct measurement of the chemical potential 
and, thus, may lead to a corroboration of some aspects here discussed. 

We proceed as follows. Section II is a brief review of general thermodynamic results 
useful for the following sections. Section III is devoted to show the relevant role played by 
the chemical potential. It is first argued that one can define it in an unambiguous way. 
Then, using its main definition, as being the negative of the change of entropy as a function 
of varying particle density at constant energy, we state the main hypothesis that negative 
chemical potentials correspond to normal phases while positive ones to quantum phases, 
zero value being the onset of the transition. A quantum phase is identified as that one in 
which there appears a condensate, as defined above, that does not contribute to the entropy. 
These phases, for interacting fluids, show the phenomenon of superfluidity. In this same 
section we discuss ideal Bose and Fermi gases and show that, besides the common ideal 
BEC, there clearly appears a phenomenon that may be called Fermi-Dirac Condensation. 
In Section IV, based mainly in experimental phase diagrams of ^He and ^He we attempt 
to build a qualitative "generic" phase diagram that obeys all thermodynamic requirements. 
The main result is the phase diagram entropy density versus particle density, where one can 
plot isoenergetic curves, exemplifying our concerns. We advance the fact that the transition 
gas to superfluid in Helium and alkaline vapours is remarkably closer to that of ideal Fermi- 
Dirac condensation. Also, by continuity, we find that there should be "supersolid" phases, 
in the sense that their chemical potential is positive, although we do not claim them to be 
the same as those supersolid phases currently discussed [251 EE]. As mentioned. Section V 
is devoted to discuss the results of recent experiments in ultracold gases, in particular, we 
argue that the density profiles measured open the door for a full determination of the phase 
diagram, for both uniform and trapped regions. We conclude with some remarks. 
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II. BRIEF REVIEW OF THERMODYNAMICS 



For a pure (non-magnetic and non-polarizable) substance and confined by a vessel of 
rigid walls of volume V, all thermodynamics may be obtained from the knowledge of the 
entropy ^S" as a function of the internal energy E, the number of atoms or molecules N and 
the volume V, namely, from the relationship S = S{E, V, N). Since S is extensive, one can 
write 

S=Vs{e,n) (1) 

where s — S/V, e — E/V and n — N/V are respectively, the entropy, energy and number 
of particles per unit of volume. Thermodynamics is obtained from the following expressions 
that summarize the First Law, 

ds — fide — adn, (2) 
where the inverse temperature p — l/ksT and chemical potential a — /i/kBT are, 

with ks Boltzmann constant. We note that a is a dimensionless variable that we shall refer 
to both this and /i as the chemical potential. The pressure p of the system is obtained via 
Euler expression, that is to say, 

— s — f3e + an. (4) 

Thermodynamics assumes that s, f3, a and p are single valued functions of e and n. The 
Second Law states that s — s{e, n) is a concave function of its variables (e, n), namely. 



— ^ < and — ^ 
oe on 
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.<0 and ^-^<0. (5) 

These conditions assert the stability of the thermodynamic states, the first inequality implies 
the positivity of the specific heat at constant volume Cy while the second one the positivity 
of the isothermal compressibility, 

C^-,3=(|)_^>0 and '--^(|)^>0 (6) 

The second expression is equivalent to the more common one kt — {'\-/n){dn/ dp)^. 

The Third Law asserts that Absolute Zero, ^ — > oo, is unattainable and, if we consider 
systems with energy spectra unbounded from above, one concludes that the temperature is 



always positive, {3 > 0[27]. The Laws are silent regarding the sign of the chemical potential a. 
We shall argue that there is a profound meaning in the sign of this quantity, typically being 
negative, while its becoming positive signals the appearance of "quantum" fluid phases. But 
before getting into that aspect, let us review briefly the characteristics of first and second 
order phase transitions. 

Phase transitions only occur in the thermodynamic limit, where S, E, V, and N become 
infinitely large, while the densities s, e, and n remain finite. At a first order phase transition 
the (originally) extensive quantities s, e and n become discontinuous while the intensive 
ones, P, a and p remain continuous. A phase transition of the second order occurs with a 
continuous change of all variables, s, e, n, /3, a and p, however, the state is strictly unstable 
in the sense that the stability conditions of the Second Law, Eqs.(|5|, are no longer negative 
but become zero. That is, — > and — t- 0. Or, equivalently, c^, — > oo and kt ^ oo 
at the transition. By well known relationships of statistical physics [28] . the former indicates 
a divergence of the density energy fluctuations, while the latter a divergence of the particle 
density fluctuations. Thus, in the strict thermodynamic limit, a system cannot be stabilized 
at a second order phase transition. For our purposes below, a second order phase transition 
implies that both the pressure p and the chemical potential a, as functions of n for a given 
value of the temperature (3, become flat at the transition. 

With the previous information plus empirical data regarding the overall characteristics 
of the phase diagram of the equation of state p = p{n, (3) one can construct the structure of 
the fundamental relationship s = s{e,n). We stress out the important fact that the equation 
of state p = p{n, (3) is not a fundamental relationship in the sense that not all the thermo- 
dynamics can be obtained from its knowledge. One needs an additional function, such as 
a = a{n,P). However, s = s{e,n) is fundamental since it does contain all thermodynamics. 

III. THE CHEMICAL POTENTIAL 

A. The non-ambiguity of the chemical potential 

As mentioned above, the sign of the chemical potential a is not restricted by the Laws 
of Thermodynamics. Moreover, there is a widespread belief that the chemical potential is 
defined up to an arbitrary additive constant. Thus, enquiring about its sign and the meaning 
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of it, may appear as a futile exercise. We now argue that such a behef is not correct and 
estabhsh that the sign of the chemical potential is of fundamental relevance. 

An alternative determination of the chemical potential to that of Eq. ([s]) is the following, 

This identity follows from the uniqueness of s as a function of e and n. The argument of the 
indeterminacy of a is that since the system is non-relativistic, the energy E can be defined 
up to an arbitrary constant C thus shifting its origin arbitrarily, E' = E + C, and because 
of Eq.([7]), yielding a chemical potential relative to such an arbitray origin of the energy. 
Let us review this carefully. First, since the thermodynamic limit must be imposed, clearly, 
an arbitrary constant C would drop out since e' = E' /V + C/V — > E/V = e. Thus, the 
arbitrary constant should be extensive, that is, C = Nc, with c indeed an arbitrary constant. 
This would certainly imply /3a' = Pa + c. We now argue that the constant c can always be 
(as it is done in practice) set equal to zero. 

Consider a realistic model of a pure fluid. It is first assumed to be non-relativistic. Thus, 
the Hamiltonian of atoms (or molecules), classical or quantum, may be written as, 

N -12 N 

H = T.^ + V^ntirl,r2, . . . , r-W) + E ^e..(r,), (8) 
i=i j=i 

where Vexti'f^) is the common confining potential, and here we are concerned with rigid- walls 
potentials only. The interatomic potential is V^„t(r*i, r2, . . . , r^) and can contain two-, three-, 
etc., body interactions. This potential can be defined, again, up to an arbitrary constant. 
However, for the energy E to be extensive, it should be true that in the limit of "infinite 
dilution", namely, in the limit \fi — r^l — oo for all pairs of particles i and j in the system, 
the potential should take the form Vint{ri,r2, ■ ■ ■ ,'nv) — ^ Nc. In all textbooks is assumed 
that the arbitrary constant is zero, c = 0. That is, one always tacitly assume that in such a 
limit the system behaves as an ideal gas (quantum or classical). We shall take this point of 
view here. This means that the energy is uniquely defined and so does the chemical potential 
a. With regard to this common convention, we can establish the meaning of the sign of a. 

B. The sign of the chemical potential 

Recall the definition of the chemical potential, Eq.(|3]). It gives the negative of the change 
of entropy with respect to a change of particles, for a given value of the energy. Thus, if 



negative (positive), it tells us that the entropy increases (decreases) with increasing number 
of particles, for a fixed value of the energy. We now resort to the well-known recipe from 
statistical physics pS] to calculate the entropy of an open system in interaction with a heat 
bath, in terms of its (average) energy E, number of particles N and volume V, 

S{E,N,V) = kBlnAr{E,N,V), (9) 

where Ar{E,N,V) is the number of states the system has access when in thermal equi- 
librium with its environment. This number of states is defined within an interval AE, the 
energy fluctuation of the thermodynamic state, namely, 

Ar=§A£, (10) 

where dT/dE is the density of states. For a fixed value of the energy E and the volume V, 
the statistical weight Ar(A^) is a function of the number of particles. 

To elucidate the dependence of s on n, let us consider a gas at large enough temperatures 
such that the system may well be approximated as classical and even close to an ideal gas. 
Let us keep the energy e fixed at that corresponding value. Then, it should be clear, and 
textbook experience corroborates it, that if the number of particles is increased the number 
of available states does too increase. That is, there is an increase of available states by 
enlarging the phase space with more particles. The chemical potential a is thus negative. 
One finds, however, that the energy per particle e = e/n obviously decreases and, at the 
same time, not completely obvious but corroborated by simple cases and verified below, the 
temperature T also decreases. It is important to realize that as the number of particles is 
increased, the structure of the states and of the energy spectrum do change, that is, the 
system is described by a different Hamiltonian, since this depends on A^. Consider now that 
the energy is still kept constant but the number of particles keeps increasing even further. 
Two situations may occur: 

i) The system reaches a limiting state of finite temperature. That is, the system reaches 
a state where if the number of particles is increased even more, the energy can no longer be 
kept constant. That is, the system must suffer a first order phase transition, the temperature, 
pressure and chemical potential changing continuously, but the energy and entropy becoming 
discontinuous, typically increasing. In brief, the isoenergetic curves eventually meet a state 
that borders an unstable region, where phase separation occurs, but the important point 
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is that the chemical potential is always negative a < since, in this case, entropy is an 
increasing function of n . We call "normal" to these states. 

ii) The system reaches the limiting state of zero temperature. This case may occur because 
the energy spectrum of the many-body system is bounded from below by the ground state 
energy Eq = Eq{N, V). Thus, for a fixed value of the energy e, there exists a maximum value 
of the density n^ax at which the energy of the ground state of that density equals the energy 
e, namely, eo{nmax) — e, where cq = Eq/V . At this stage the system is in its corresponding 
ground state, the entropy becomes zero s(e, Umax) — and so does the temperature T — Q. 
Since at low enough density the entropy increases with increasing density but eventually 
reaches s = for rimax, it follows that the entropy must have reached a maximum value at 
some critical density n*, above which the entropy begins to decrease. This in turn implies 
that the chemical potential is negative for small enough density, becoming zero at n* < Umax 
and turning into a positive quantitive above it. We shall argue further below that this change 
signals the appearance of, or transition to, a quantum fluid or a quantum solid with proper 
macroscopic properties. Let us advance a hypothesis of the nature of that state that explains 
the decrease of the entropy at the point a = 0. 

C. The nucleation of a condensate phase 

As mentioned above, if for a fixed value of the energy the system reaches a limiting state 
at zero temperature, which in turn is the ground state for that density, it must be true 
that in the vicinity of that state the chemical potential must be positive indicating that 
the entropy decreases with increasing density. This implies that at some density n* the 
chemical potential becomes zero Q;(e, n*) — 0. One can conclude that as this point is crossed 
by increasing the density, a macroscopic phase nucleates such that it does not contribute 
to the entropy. To have such a property, this phase must be a single many-body quantum 
state. As we shall analyze below, this coincides with the appearance of supcrfluid phases 
in He and alkaline ultracold vapors. For a Bose ideal gas, this is also the phenomenon of 
Bose- Einstein Condensation (BEC). What we further claim here, as the main hypothesis 
of this study, is that this phase must be very close to the many-body ground state of the 
fraction of the particles that conforms such a phase. Let Uc be the fraction of particles (per 
unit volume) in the zero-entropy phase and ut — n — Uc the remaining ones. Thus, we claim, 
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the zero-entropy phase is the many-body ground state of the Uc particles with energy eo(nc), 
i.e. e(n) = exinT) + eo(?T.c). As the number of particles is increased, increases up to the 
point where n^- = and Uc = Umax, ssid the whole system is in its corresponding ground 
state. The entropy and the temperature are zero. We immediately point out that this is not 
the usual Bose-Einstein condensation scenario. In such a case the condensation occurs at a 
sm^'/e-particle state. Here, we insist in the fact that the condensate is the many-body ground 
state of the condensate particles. Of course, the many-body ground state of ric particles 
per unit volume in the ideal Bose gas is the same as particles per unit volume in the 
ground state of one particle. We discuss now that the most common case of condensation of 
a many-body ground state is already present in the ideal Fermi gas and that this is actually 
closer to the scenario of real superfluids. 

D. Ideal Bose-Einstein and Fermi-Dirac Condensation 

It is common knowledge that Bose-Einstein condensation in ideal gases is a phenomenon 
characteristic of Bose statistics only. Here we argue that as a matter of fact, the phenomenon 
of condensation, understood as the development of a quantum phase that decreases the 
entropy of the system and characterized by a positive chemical potential, is also present in 
ideal Fermi-Dirac gases. In Figs. [T] and [2] we show the entropy s as a function of the density 
n for fixed values of the energy e, for ideal gases[29]. Fig. [l]is for Fermi and Fig. [2] for Bose 
statistics. Both curves begin, at very low density, with chemical potential negative. At some 
density n* both curves reach a maximum and the chemical potential a becomes zero. For 
bosons this is the onset of the well known phenomenon of Bose-Einstein condensation, and 
beyond this point the chemical potential remains zero. The interpretation is that the single 
particle ground state begins to be macroscopically occupied. That is, at fixed energy, for any 
given value of the density n > n*, there is a fraction Uc = n — n* that occupies the single- 
particle ground state. We note, however, that such a state is also the many body ground 
state of the particles. Moreover, the energy of the system equals that of the critical 
density e = e(n*) since the condensate has zero energy. For fermions, the zero chemical 
potential point is essentially ignored in most textbook discussions. We argue now that this 
also marks the onset of a phenomenon that may also be called Fermi-Dirac condensation. 
The condensation is not the macroscopic occupation of a single-particle state but rather of 
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FIG. 1. Entropy s versus density n at a fixed value of the energy e = 0.5, 1.0 and 2.0 for an ideal 
Fermi gas, in units Ti = m = ks = ^- The slope of the curve is the negative of a = /3//. 

a many-body ground state, that of the condensed fraction. 

From Fig. [T]we see that, at fixed energy, the entropy reaches a maximum at a density 
value n* where the chemical potential becomes zero. Beyond this point the chemical potential 
becomes positive indicating a decrease of the entropy. What does occur in the gas? A 
macroscopic occupation of a single-particle state is forbidden by Fermi statistics. However, 
we argue, there appears a "condensate" fraction which is in a single many body state in 
order to give a zero contribution to the entropy. This condensate grows a n is increased thus 
decreasing the entropy. This behavior continues up to a maximum value Umax where the 
entropy is zero and where all the gas is in its corresponding many body ground state. The 
latter state has an energy eoirimax)- Note that at this point a — )■ oo, /3 — )■ oo but the ratio 
fi = a//3 equals the so-called Fermi energy ep. We argue then, that for values between = 
and /i = ei?, a fraction must occupy a single many-body state that does not contribute to 
the entropy. By continuity in the neighborhood of zero temperature, that is, at Umax where 
all the gas is indeed in its many body ground state, the condensate must also be in its own 
many body ground state. To find this condensate fraction we appeal to Fig. [3] where we 
plot /X vs n for different values of the temperature. The thick line denotes the value of the 
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FIG. 2. Entropy s versus density n at a fixed value of the energy e = 0.5, 1.0 and 2.0 for an ideal 
Bose gas, in units h = m = kB = 1- The slope of the curve is the negative of a = /3/i. The dotted 
lines correspond to Bose- Einstein condensation states. 

chemical potential at zero temperature fio{n) = fi(n, T = 0), i.e. fio{n) = epin). Take a fixed 
value of n, say the vertical dotted line in Figjsj Let Tc(n) be the temperature at which the 
chemical potential is zero. For temperatures T < Tc{n), the chemical potential is positive, 
/i > 0, and the total fixed density equals n = tit + n^^n), where no(/i) is the value of the 
density obeying fJ.o{no) = /i. This is indicated by the horizontal dotted line in Fig. [3] The 
thermal density ut is obtained from the difference of n — nQ{^). For temperatures T > Tc{n), 
clearly no(/u) = 0, all the gas is in a thermal state and the chemical potential is negative, 
/i < 0. Figure |4] is a plot oiriQ/n vs T/Tc{n) obtained from Fig. [3] for any value of n, that is, 
it is a universal curve. In analogy to Bose-Einstein condensation, this phenomenon may be 
called Fermi-Dirac condensation for particles obeying Fermi statistics. We repeat once more, 
the condensate is the fraction of the density that is in its own many-body ground state. As 
explained above this state coincides with the one-particle ground state condensate for ideal 
Bose gases. In the next sections we shall present a generic phase diagram for real interacting 
quantum fluids that present the phenomenon of superfluidity. We will see that such a state 
behaves more like the ideal Fermi-Dirac condensation than the ideal Bose-Einstein one. 
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FIG. 3. Isotherms (T = const) fivsn for an ideal Fermi-Dirac gas. The thick curve is T = (Fermi 
energy). The vertical dashed line indicates an arbitrary density n. This line crosses all isotherms 
and defines the critical one Tc(n) at ^ = 0. For temperatures T < Tdn) the density separates 
into a condensate fraction no and a thermal one n-r- at the temperature T defined by the crossing 
of the vertical and horizontal dashed lines, the condensate fraction no(r)is the intersection of the 
horizontal dashed line with the T = curve, while the thermal part is nxiT) = n — no(r). Clearly, 
no becomes zero for T > Tc{n) and no = n at T = 0. 



IV. GENERIC PHASE DIAGRAMS OF SUBSTANCES THAT SHOW QUAN- 
TUM PHASES 

As mentioned in Section II, the relationship s = s{e,n) is a fundamental one in the sense 
that it contains all the thermodynamic information of the system. To construct it we may 
follow the recipes of statistical mechanics [28] and/or resort to empirical evidence. The latter 
is typically obtained with the measurement of the equation of state p = p{n,T). In this 
section, based on the empirical equations of state of ^He and ^He[Hl EHl EI], and using the 
general properties of the equilibrium states as describe in Section II, we construct generic 
phase diagrams of s = s{e,n). In the next section, we argue that many of the characteristics 
of those phase diagrams should remain valid for bosons or fermions systems whether confined 
or not. 
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FIG. 4. Fermi-Dirac Condensation. Plot of no/n vs T/T^ where nQ{n,T) is the fraction of particles 
that occupy its own many body ground state. Tc{n) is the critical temperature for density n 
obtained when the chemical potential becomes zero, see Fig. [Sj 




FIG. 5. Qualitative phase diagrams of a pure substance showing a superfluid phase. The areas in 
gray color are forbidden unstable regions. 
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Fig. |5]\ shows a phase diagram p vs T that quahtatively resembles that of ^He and 
^He[Hl EHl EI]- Four phases are present, gas, hquid, superfluid and sohd. The hquid - sohd, 
hquid - gas and superfluid - solid transitions are believed to be first order phase transitions, 
while liquid - superfluid and gas - superfluid second order (as well as the isolated critical 
point in the normal liquid - gas transition). We consider all first order transitions being 
discontinuous in the particle, energy and entropy densities. The transition line liquid to 
superfluid is supposed to have positive slope. This is true in ^He while it is negative in 
^He .A transition line with negative slope is considered "anomalous" since some properties 
do not follow "expected" behavior near the transition. However, there is no fundamental 
requirement for the sign of coexistence lines and we have assumed all of them to be positive 
for simplicity. We mention that even at zero magnetic field the ^He phase diagram shows 
two superfluid phases (so called A and -B)[3l [20] due to the particular magnetic properties 
of ^He. Again, for simplicity, we consider our hypothetical system to show the simplest of 
phase diagram, including nevertheless, a superfluid phase 

Based on Fig. [5]4. and on general thermodynamic requirements one can construct the 
qualitative phase diagram p vs n. This is shown in Fig. [5^. In Figs. [6] and [7] we show a 
couple of isotherms obeying the general requirements of Section II, indicating how the phase 
diagram p vs n of Fig. |5]A was constructed. On the one hand, we appeal to the positivity 
of the isothermal compressibility, Eq.(|5|, to obtain p as an increasing function of n for T 
constant. Second, we observe that in first order phase transitions n is discontinuous with p 
continuous, while p shows a flat slope at the density n of second order transitions. 

The relevance of the phase diagram p vs n shown in Fig. [5] is that its allows us to qual- 
itatively construct the phase diagram s vs n. With this, and following the thermodynamic 
laws, one can then built s vs e. These diagrams are shown in Fig. |8j These two diagrams 
embody all the thermodynamics of the system. For the purposes of the present article, we 
show in Fig. [9] different isoenergetic lines s vs n whose meaning is discussed below. The 
discussion of the isodensity curves s vs e is not presented in this article since they can be 
very complicated and do not add more to the main point of this section. 

Returning to Fig. [9} we show several isoenergetic lines s vs n. By the stability of the 
states, see Eq.([5]), their curvature must be negative and its slope is a = /3yU. Since /3 is 
always positive, this implies that the chemical potential has the same sign as a. In Fig. [9] 
we have indicated the regions of the phase diagram where fi is either positive or negative. 
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FIG. 6. Isotherm showing a gas-hquid first order, hquid-superfluid second order and superfluid- 
soUd first order phase transitions. By stabihty, see Eq.([5]), the slope of the isotherm p ys n is 
always positive, except at a second order phase transition where it becomes flat (marked with a 
dot) showing the critical divergence of the compressibility. The areas in gray color are forbidden 
unstable regions. 



We note that a fortiori there are regions where the chemical potential is negative and others 
where is positive. That it must be negative at some points is a consequence that in the 
appropriate limit the normal gas behaves as an ideal gas, while the positive sign region 
must exist because very near to T = the system is essentially in its ground state and the 
entropy must be reduced as n increases for e fixed. We note the striking similarity of the 
isoenergetic curves in the gas to superfluid region of Fig. [9] with that of the ideal Fermi gas 
in Fig. [Tj Our claim that this region marks the existence of a condensate fraction, being a 
many body ground state as explained above, is supported by the fact that the region /i > 
coincides with the superfluid macroscopic phase state. That is, it is well accepted that the 
superfluid fraction of the superfluid state at T 7^ does not contribute to the entropy. It 
appears natural to identify the condensate as the superfluid fraction. 

We would like to add that due to the continuity of the chemical potential (and temper- 
ature and pressure) at a flrst order phase transition, then, there must be solid phases with 
chemical potential negative and others with positive. The former are the "normal" solids, 
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FIG. 7. Isotherm showing a gas-superfluid second order (marked with a dot) and a superfluid-sohd 
first order phase transitions. The areas in gray color are forbidden unstable regions. See caption 
of Fig. 6. 

while the latter indicate the existence of a "quantum" solid or "supersolid" . The transition 
line between these two is within the solid and it should be a line of zero chemical potential. 
This is also indicated in the figure. We are not claiming that the "supersolid" phase is 
necessarily the same as that recently discussed in the literature [25t [26]. Further analysis of 
this phase is beyond the scope of this article. 

To conclude this section, we point out that if we had started with a phase diagram 
p = p{n, (3) of a pure "normal" substance, that is, one that does not show superfluid phases, 
one obtains that the chemical potential is always negative. 

V. MANY BODY THEORIES AT T = AND RECENT EXPERIMENTS IN UL- 
TRACOLD CONFINED GASES 

As mentioned in the Introduction, many-body theories have predicted quite satisfactorily 
the ground state of interacting gases, thus describing the superfluid state at zero tempera- 
ture. On the other hand, the details of the phase transitions to superfluid states is still a 
challenge to theoretical efforts. From experimental evidence is believed that transitions to 
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FIG. 8. Phase diagrams of the fundamental relationship s = s{e,n). The areas in gray color are 
forbidden unstable regions. 

superfluid from gas and liquid phases are second order. However, even at the level of mean 
field theories the description of those transitions is still a matter to be settled [321 [33]. In 
addition to the predictions of the celebrated theories mentioned above, we believe that the 
present experiments with ultracold gases show the opportunity to prove, if not inadvertendly 
already been done, that the phase diagrams presented above are the correct ones. In partic- 
ular, the ultracold experiments have probed the region of very small temperatures and very 
small pressures, where the transition gas to superfluid occurs. In this section we address the 
phenomena observed in trapped inhomogenous ultracold gases, to argue that an analysis of 
experimental data can be made that should lead directly to the appropriate phase diagrams 
yielding further information on the order of the transition. 
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FIG. 9. Phase diagram s vs n showing isoenergetic curves indicating the regions of normal fi < 
and quantum fi > behavior. Stabihty requires those curves to have negative curvature. Its slope 
is the negative of a = We identify the normal to quantum transition as being characterized 
by /X = 0. 

For our purposes we recall first the predictions of the fundamental theories regarding 
the behavior of the chemical potential as a function of the density. At zero temperature, 
both for fermions and bosons, the chemical potential is a positive increasing function of the 
density. That is, the chemical potential predicted by Bogolubov and BCS is very similar to 
that of an ideal Fermi gas (the thick line in Fig. [sj. We recall that Bogolubov prediction is 
/i(T = 0) ~ n, while for BCS is essentially the ideal Fermi energy fi{T = 0) ~ n"^^'^. This 
result, in addition to the fact that superfluidity is present, seems to be in accord with the 
ideas presented in the previous sections. That is to say, at some point the chemical potential 
must have changed sign indicating the appearance of a condensate-superfluid phase. 

A very important piece of information is the way in which the transition is made. If second 
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order, the chemical potential should be flat at the transition, identified as the point where 
the chemical potential changes sign. Let us consider the isothermal functions of chemical 



potential versus density. In Fig. 10 we show an actual curve jj, vs n for a temperature 



T = 1.0 for an ideal Fermi gas, while in Fig. 11 we show a hypothetically expected isotherm 



for an interacting gas, be it Fermi or Bose. The main difference between the latter and 
the former is that the slope of /i vs n at the transition /i = is zero, as dictated by the 
fact that it should be a second order phase transition. The tail for /i < should become 
that of an ideal classical gas, while the corresponding one for > should follow that of 
Bogolubov or BCS theories, as the temperature nears zero. The difference of those curves 
being quantitative but certainly not qualitative. Experiments with ultracold gases have the 
peculiarity that are performed in confining potentials, typically harmonic, that give raise 
to inhomogenous gases. It is well established that the so-called local density approximation 
(LDA) yields accurately the density profile measured in these experiments [3H439] . This 
approximation indicates that the density profile, for an isotropic harmonic potential, may 
be calculated as 

p(r) = n{fi — -mw^r^, T). (11) 
That is, the equation of state n = n{fi,T) for the homogenous system is the inhomogenous 



density profile rescaled by the confining potential. Fig. 12 is the profile for an isotherm of 



the hypothetical interacting system of Fig. 11, for a temperature and chemical potential of 



the confined system below condensation. A very important prediction of the present study 
is thus found: the point where the chemical potential is zero in the homogenous system 
translates to the point where the profile changes behavior. Call this spatial point Tc, then 
the chemical potential of the confined inhomogenous gas is yU = ^mcu^r^. That is, density 
profiles are enough to construct not only the phase diagram of the homogeneous systems but 
also of the inhomogenous ones. The latter follow from the fact that the number of particles 
in the confined gas is 

N = f p{r)d\. (12) 



This yields = N{T, fi, u~'^) which is the equation of state of the confined gas as a function 
of T, fi and the "generalized volume" o;"^ [351 [36]; see also Ref. |10] where an effort has been 
made to experimentally construct the phase diagram of a trapped ^^Rb gas. Although finite 
size effects may obscure the behavior of the density profiles at the region where they typically 
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FIG. 10. Isotherm n vs /i for T = 1.0, of an ideal Fermi gas. Units h = m = ks = 1- 

change curvature, systematic measurements with different number of particles should help 
to elucidate this relevant point. 

We recall here that building phase diagrams from density profiles has already been 
performed [371439] . however, the identification of zero chemical potential as the transition 
points has not been used. 

VI. FINAL REMARKS 

In his celebrated book [19] on superfluidity, Fritz London makes a description of the na- 
ture of the superfiuid phase that essentially has survived up to these days, and that it is 
also the main point of the present article although with a twist. London was impressed 
by the fact that the superfiuid fraction does not contribute to the entropy of the fluid. He 
considers the fraction to be a single macroscopic quantum state and he even says that su- 
perfiuid flow allows to "tap ground state" from the fluid[19]. The further but inaccurate 
London's assumption is that the superfiuid state may be considered to be similar to an ideal 
Bose- Einstein condensate. It is well recorded that Landau [10] argued that it could not be 
so since an ideal gas cannot be superfiuid. Few years later, Bogolubov[TT] provided the first 
satisfactory model of superfluidity at T = 0, which needed the presence of intermolecular 
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FIG. 11. Hypothetical isoterm n vs fi of an intereracting quantum gas. As a matter of fact, for 
/i < we used an ideal Bose gas and a curve ^ = C{n — Uc)^ with 7 > 1, for fj, > where C is a 
constant and ric is the "critical" value of the density, namely, its value at /u = 0. The hypothesis 
is that the transition should be a second order phase transition and, therefore, the slope of ;U vs n 
should be zero at fi = 0. 

interactions in an essential way. Nevertheless, since then, many of theoretical efforts and 
explanations have been prone to describe superfluidity as a kind of Bose-Einstein conden- 
sation in a single particle state, although certainly including interactions. It may seem, as 
London himself believed, that Bose statistics are essential. However, we also well know that 
this is not necessarily true. BCS-like theories have been successfully applied f3j to explain 
superfluidity in ^He, a fermion. Thus, the existence of quantum fluids phases that may 
exhibit superfluidity is not truly a matter of statistics, it is a phenomenon due solely to the 
quantum nature of the fluids and to the fact that the samples are macroscopic. Of course, 
the details may depend on the particular fluid, and thus certain properties may depend on 
being bosons or fermions. And it is in this regard that the main thesis of this article rests: 
the phenomenon of condensation is the sudden occupation of certain number of atoms of 
their many-body ground state, hence not contributing to the entropy, and reaching the full 
ground state of the whole sample at zero temperature. Clearly, this scenario includes ideal 
Bose-Einstein condensation, but it does appear that ideal Fermi-Dirac condensation is closer 
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FIG. 12. Hypothetical density profile p{r) vs r of an interacting quantum gas confined by a 
harmonic trap, for a temperature below condensation T < Tc- This profile was obtained using the 



local density approximation, see Eq.( 11 ), of the isotherm /i vs n in Fig. 11 The frequency is w = 1 
in units of fi = m = ks = 1- Typical experimental density profiles look like this, with a thermal 
gaussian tail and a peak indicating the presence of the condensate. 



to the interacting case, be it fermions or bosons. One of the suggestions of this article is 
that such an observation may be useful in the search of a complete and satisfactory theory 
of the thermodynamics of quantum fluids. 

ACKNOWLEDGMENTS 

We acknowledge support from grant PAPIIT-UNAM IN116110. I thank Dulce Aguilar 
for technical assistance in the preparation of the figures of this article. 



[1] P. Kapitza, Nature 141, 74 (1938). 

[2] J.F. Allen and A.D. Misener, Nature 141, 75 (1938). 

[3] D. D. Osheroff, R. C. Richardson, and D. M. Lee, Phys Rev Lett 28, 885 (1972) 

22 



[4] For a recent account of the discovery of superfluidity in ^He, see S. Balibar, J. Low Temp. 
Phys. 146, 441 (2007). 

[5] M. H. Anderson, J. R. Ensher, W. R. Matthews, C. E. Wieman, and E. A. Cornell, Science 
269, 198 (1995). 

[6] K. B. Davis, M. -O. Mewes, M. R. Andrews, N. J. van Druten, D. S. Durfee, D. M. Kurn, and 
W. Ketterle, Phys. Rev. Lett. 75, 3969 (1995). 

[7] C. C. Bradley, C. A. Sackett, J. J. ToUett, and R. G. Hulet, Phys. Rev. Lett. 75, 1687 (1995). 

[8] F. London, Nature 141, 643 (1938). 

[9] L. Tisza, Nature 141, 913 (1938). 
[10] L. D. Landau, Phys. Rev. 60, 356 (1941). 
[11] N. N. Bogolubov, J. Phys. USSR 11, 23 (1947). 

[12] J. Bardeen, L. N. Cooper, and J. R. Schriefer, Phys Rev 108, 1175 (1957). 
[13] A. J. Leggett, Phys. Rev. Lett. 29, 1227 (1972). 

[14] M. R. Matthews, B. R Anderson, P. C. Haljan, D. S. Hall, C. E. Wieman, and E. A. Cornell, 

Phys. Rev. Lett. 83, 2498 (1999). 
[15] K. W. Madison, F. Chevy, W. Wohlleben, and J. Dalibard, Phys. Rev. Lett. 84, 806 (2000). 
[16] J.R. Abo-Shaeer, C. Raman, J.M. Vogels, and W. Ketterle, Science 292, 476 (2001). 
[17] For a recent review of vortices in BEC, see A. L. Fetter, Rev. Mod. Phys. 81, 689 (2009). 
[18] M. W. Zwierlein, Abo-Shaeer, J. R., Schirotzek, A., Schunck, C. H. and Ketterle, W. Vortices 

and superfluidity in a strongly interacting Fermi gas. Nature 435, 1047 (2005). 
[19] F. London, Superfluids, John Wiley &; Sons, New York (1954). 
[20] E. R. Grilly and R. L. Mills, Ann. Phys. 8, 1 (1959). 

[21] B.D. Holian, W.D. Gwinn, A.C. Luntz, and B.J. Alder, J. Chem. Phys. 59, 5444 (1973). 
[22] J.K. Hoffer, W. R. Gardner, C. G. Waterfield and N. E. Phillips, J. Low. Temp. Phys. 23, 63 
(1976). 

[23] A. Driessen, E. van der Poll, and LF. Silvera, Phys. Rev. B 33, 3269 (1986). 
[24] Y.H. Huang, G.B. Chen, B.H. Lai, and S.Q. Wang, Cryogenics45, 687 (2005). 
[25] M. H. W. Chan, Science 319, 1207 (2008). 
[26] S. Balibar, Nature 464, 176 (2010). 

[27] If the energy spectrum is also bounded from above, there are negative temperatures, see N. 
F. Ramsey, Phys. Rev. 103, 20 (1956). 

23 



[28] L. D. Landau and E. M. Liftshitz, Statistical Physics I , Pergamon Press, London (1980). 
[29] All the calculations regarding ideal gases follow straightforwardly from the expressions for the 

grand potential, as found in any textbook in Statistical Physics, see for instance Ref. [28] 
[30] A. M. Guenault, Basic superfluids (Master Series in Physics and Astronomy, CRC Press, 

London (2002). 



[31] The web page http://ltl.tkk.fi/research/theory/he3. html of the Research on helium of the 
LTL/Helsinky University of Technology shows phase diagrams of '^He under several condi- 
tions. 

[32] L. Olivares-Quiroz and V. Romero-Rochm, J. Phys. B: At. Mol. Opt. Phys. 43, 205302 (2010). 
[33] L. Olivares-Quiroz and V. Romero-Rochm, J. Low Temp. Phys. 164, 23 (2011). 
[34] F. Dalfovo, S. Giorgini, L. P. Pitaevskii, and S. Stringari, Rev. Mod. Phys. 71, 463 (1999). 
[35] V. Romero-Rochm, Phys. Rev. Lett. 94, 130601 (2005) 

[36] N. Sandoval-Figueroa V. and Romero-Rochm, Phys. Rev. E 78, 061129 (2008). 
[37] M. Horikoshi, S. Nakajima, M. Ueda, T. and Mukaiyama, Science 327, 442 (2010) 
[38] T.-L Ho and Q. Zhou, Nature Phys. 6 131 (2010). 

[39] S. Nascimbene, N. Navon, K- J. Jiang, F. Chevy F, and C Salomon, Nature 463 1057 (2010). 
[40] V. Romero-Rochm, R. F. Shiozaki, M. Caracanhas, E. A. L. Henn, K. M. F. Magalhaes, G. 

Roati, and V. S. Bagnato, BEC phase diagram of a ^"^ Rb trapped gas in terms of macroscopic 

thermodynamic parameters, Cond-mat arXiv:1102.4642| 



24 



